Interspecies singlet pairing in a mixture of two spin-1 Bose condensates 
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We study the ground state properties of a mixture formed by two spin-1 condensates in the 
absence of an external magnetic field. As the collisional symmetry between interspecies bosonic 
atoms is broken, the interspecies coupling interaction (/?) and interspecies singlet pairing interaction 
(7) arise. The ground state can be calculated using the angular momentum theory analytically for 
7 = 0. The full quantum approach of exact diagonalization is adopted numerically to consider the 
more general case as 7 7^ 0. We illustrate the competition between the two interspecies interactions 
and find that as singlet pairing interaction dominates (or the total spin vanishes), there are still 
different types of singlet formations which are well determined by j3. 
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I. INTRODUCTION 

Since the MIT group succeeded in trapping a ^^Na 
condensate in an optical potential [l|, the spin degrees 
of freedom are liberated, that give rise to a rich variety 
phenomena (^] such as spin domains Q , textures [3| , spin 
mixing dynamics 0-0] , and fragmentation of condensate 
[1, [13 • The properties of such a three-component spinor 
condensate were first studied using mean-field theories 
(MFT) 0, and polar and ferromagnetic spinor conden- 
sates had been implemented experimentally [ll|, US] • 1^ 
was initially predicted that the ground state of ^'^Na BEC 
(c2 > 0) is either polar (no — N) or anti-ferromagnetic 
(rii — n_i ~ N/2) in the mean-field theory. However, 
the results from many body theory Q pointed out that 
the ground state of ^^Na atoms is a spin singlet with 
properties drastically different from those of mean field 
theories (ni — no, = n^i ~ N/3). It was shown Q 
that the singlet ground state in zero field is a fragmented 
condensate with anomalously large number fiuctuations 
[Aria ^ N) and thus has fragile stability. The exact 
quantum eigenstates in the spin-2 case [13, are also 
found and their magnetic response to a weak magnetic 
field is compared with their mean field counterpart [l5l |. 

Mixtures of scalar condensates with more than one 
atomic species or state are actively studied theoretically 
[T6l - [T9| . Experimentally, the Feshbach resonance has ex- 
ploited to create a double species condensate with tun- 
able interactions and the dynamics of the superfluid and 
controllable phase separation are observed [20|-[23j. By 
adjusting the two s-wave scattering lengths ap and 02 
through the so-called optical Feshbach resonances [23 |. 
the spin exchange interaction between individual atoms 
can be precisely tuned. The theoretical studies on mix- 
tures of spinor condensates attracted much attention re- 
cently and both mean field and quantum many body the- 
ories have been applied to this novel system |25l - l29l [. A 
temporal modulation of spin exchange interaction, which 



is tunable with optical Feshbach resonance, was recently 
proposed to localize the spin mixing dynamics in a ®^Rb 
condensate [soj . 

The interspecies scattering parameters between *^Rb 
and ^'^Na are calculated resorting to the simple approach 
of the degenerate internal-state approximation (DIA) 
[31I - I33I as the low-energy atomic interactions can be 
mostly attributed to the ground-state configurations of 
the two valence electrons. The interspecies scattering 
lengths for singlet and triplet electronic states are ap- 
proximately determined already [13, 113 , given by 05 = 
109ao and ot = 70ao, where ag is the Bohr radius. The 
interspecies interactions between ^''Rb and ^'^Na atoms 
are then parametrized by three scattering lengths (see 
eq. Il]) bellow), each being a linear combination of as 
and ax weighted by the appropriate 9j coefficients for 
the total combined spins of _F = 0, 1, and 2. Within this 
approximation, it is found coincidcntly that, the param- 
eter for the interspecies singlet-pairing interaction 7, is 
equal to zero ^9| . 

However, for the more general case of an arbitrary 
mixture of spin-1 condensate, the hyperfine interaction 
between nuclear spin and electron spin gives in general 
non-neglegible interspecies pairing and DIA approxima- 
tion is not applicable [1^ Isil - fssj . In this paper, we first 
study the various quantum phases of the binary mixture 
of spin-1 condensates in the ground state ignoring the in- 
terspecies singlet pairing. Then the situation with com- 
petition between intra- and inter-species singlet pairings 
is considered. Using the full quantum approach of exact 
diagonalization, we present the detailed phases diagram 
for the 7 7^ case. 



II. THE HAMILTONIAN 

We take the intra-condensate atomic interaction in the 
form Vj{r) = (aj + /SjFj ■ Fj)S{r) with j = 1,2 for the 
two species and the inter-species interaction is described 
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In contrast to intra-condensate interactions between 
identical atoms [1^ [1^ , the collision between atoms be- 
longing to different species can occur in the total spin 
F = 1 channel, which makes the mixture more interest- 
ing. Here 
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with Cq 12 is the scattering lengths in the channels of 
total spin F 0,1, and 2, respectively, and fi = 
M1M2/ {Ml + M2) denotes the reduced mass for the 
pair of atoms, one each from the two different species 
with masses Mi and M2 respectively. ^0,1,2 is the 
corresponding projection operator with the relationship 
1 ^ P2+P1+P0 and Fi -Fa = P2-Pi~2Pq, from which 
we get 



-l3Yi-Y2+iP^)5{v) 
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with the parameters a = {g'^i'^^ +92^^' = {"dY 

5^^') )/2 and 7 = (2g<''^ - 85^'^ )/2^o projects 

an inter-species pair into spin singlet state [23 • Further- 
more, these parameters can be related to the singlet and 
triplet scattering lengths by means of a method based on 
9j coefficient [29j 
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Denote the atomic field operators for the spin state 
as for species 1 and $i for species 2, the Hamil- 
tonian for the mixture system in the second quantization 
is represented by 



H 
Hi 



Hi + H2 + -ffl2 
ft2 



(5) 



H12 = dr{a^^^,^, 



(6) 



H2 is the same as Hi with the substitution of subscript 
1 by 2 and by and O = - #0*0 + ^-1^1. 

Through the control of the trapping frequency, we can 
make the two species sufficiently overlapped and adopt 
the single spatial-mode approximation (SMA) [1, H, Isj] 
for each of the two spinor condensates with modes ^'(r) 
and $(r), i.e. 



(7) 



with fli (bi) the annihilation operator for the ferromag- 
netic (polar) atoms satisfying [0^,0^] = and ai,aj 

dij (and the same form of commutations for bi). The 
density-density interaction part is a constant. Hence we 
only focus on the spin-dependent Hamiltonian 



H 



ClPl 



2 

C12/3 



{¥\ ~ 2Ni) 



^(F? 



2N2) 



2 Fi .1^2-1-^612912, 



(8) 



with Fi = alFiijCLj (F2 = 6,]F2ij6j) defined in terms of 
the 3x3 spin-1 matrices Fi^ {Y2ij)- The operator 
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creates a singlet pair with one atom each from the two 
species, in much the same way that the operators 



it^(aS)2-2aIaLi, 
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create pairs with two atoms from the same species [3, [3]. 
The interaction coefficients are ci — / dr |5'(r)|^ , C2 = 
/ dr |$(r)|^ and C12 = / dr |*(r)|' |$(r)|'. 

In the following we perform all analysis in the single- 
mode regime in the absence of external fields. As learned 
from previous studies [s^l , SMA is shown to be exact for 
atomic interaction of the ferromagnetic type. For polar 
interaction, if the magnetization M = 0, the SMA wave 
function is still exact and becomes invalid only if M is 
large. In our mixture, we can safely apply the SMA to 
wave function of both polar and ferromagnetic atoms in 
the absence of an external magnetic field. The ground 
state for polar atoms is a fragile fragmented state with 
M = 0. Small external fields or spatial dependence drive 
the system to symmetry broken states which are better 
captured by mean field theory. Very recently the symme- 
try broken phase has been found in this binary mixture 
in the presence of a weak magnetic field [2^ . 

In our previous study [27[ we report the anomalous 
fluctuations for the numbers of atoms in the mixture of 
^■^Na (polar) and ®''Rb (ferromagnetic) condensates in 
their F = 1 manifold. DIA has been adopted to ignore 
the 7 term. The ferromagnetic ^^Rb condensate provides 
a smooth background where the quantum many body 
states are hardly affected by the fluctuation, while the 
fragile polar atoms ^'^Na are easier influenced. Especially 
in the ground state of the AA phase, the interspecies anti- 
ferromagnetic spin-exchange is large enough to polarize 
both species and a maximally entangled state is realized 
between two species with total spin F — 0. 

III. THE SITUATION OF 7 = 

In this section, we consider the simple situation of 7 = 
with Hamiltonian 
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^t? + ^Fl + ^Fi.F2. (11) 
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or in an alternative form 



H = aFl + bFj 



cF' 
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Ci2/3/4,fo = C2/32/2 - Ci2/3/4,c = 
- F2 is the total angular momentum 
A constant energy shift of CjPjNj has been 



where a = ci/?i/2 
C12/3/4, and F = Fi 
operator 

trivially eliminated. The results may serve as reference 
states for the complete ground-state phases. 

The eigenstates of the Hamiltonian ([T2|) are then sim- 
ply the common eigenstates of the commutative opera- 
tors F?, F2: and F^, denoted by 



\Fi,F2,F,m) 



|Fi,mi) 1-^2,^2) 

(13) 



where the states in uncoupled representation 

|Fi,mi) = Z;-'{Fi^f^-"'^{alf^{A^)'^'''-^''>^^\0) 

|F2,m2) - Z2"5(/.2_)F2-™2(^t)i=^2(^t)(JV2-F2)/2 |Q^, 

span a Hilbert space of dimension (iVi_2 4- l)(^i,2 + 2)/2, 
respectively [Ij]. Here C's are the Clebsch-Gordon 
coefficients, Zi^2 are the normalization constants and 
-F'i_(i^2-) is the lowering operator for mi(m2). Minimiz- 
ing the corresponding eigenenergies of the Hamiltonian 



E = aFi{Fi + 1) + 6F2(F2 + 1) + cF{F + 1) 



(14) 



we can get the ground state energy determined by dif- 
ferent parameters ci/3i, C2/32 and ci2/3. In Fig. 1, us- 
ing the full quantum approach of exact diagonalization 

we calculate ground state order parameters (^^1'^ ■ /f''^ 



and ^Fi • F2y to illustrate the different phases as ci2/3 
changes. The results are presented for three typical cases 
(a) ci/3i = -1, C2/32 - -2, (b) ci/3i = 1, C2/32 = 2 and (c) 
ciPi = —1,02^2 ~ 2 with equal atomic numbers in two 
species, i.e. Ni — N2 — N = 100. We see that the results 
agree fairly well with the simulated annealing approach 
in the mean field theory (2 Si] except for some small devia- 
tions originated from pure quantum effect. For instance, 
the maximum and minimum values are N(N+1) — 10100 

and for (fA and (f^V but they are = 10000 and 
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-10100 for (Fi • F 
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FIG. 1: (Color online) The dependence of ground-state order 
parameters on ci2/3 at fixed values of (a) ci/3i — —1, C2P2 = 
-2; (b) ci/3i = 1, C2/32 = 2; (c) ci/3i = -1, C2P2 = 2; and 
C127 = (all in units of |ci/3i|). Black solid lines, red dashed 
lines and blue dot-dashed lines denote respectively the order 



parameters (Fi 



F^ ) and ( Fi ■ F 



with m = 0, ±1, ... ± F. The extreme states 

\F,,F2,F, F) - Cl^^}}^^j, \N, N) \N, N) (16) 

can be simply described as Z'^/'^{a[)^ {b\)^ |0) . The AA 
phase is a singlet |A^, TV, 0, 0) with all states obey the 
condition mi + 1x12 = 0. All channels of total spin zero 
have to be taken into account and we have 



|7V,7V,0,0) = 



N 



mi = -N 



(17) 



For the mixture of two ferromagnetic condensates, 
there are generally two phases FF and AA separated by 
the critical point as shown in Fig. la. The FF phase 
is described by a set of degenerate states generated by 
repeatedly applying the lowering operators (Fi_ -I- F2-) 
on the extreme states for {F — m) times 

|Fi,F2, i^,m) = {Fi^+F2-f-'''\FuF2,F,F) , (15) 



B. The case ci/3i > 0, £2^2 > 

The mixture of two polar condensates allows for five 
distinct phases FF, MM_, PP, MM+, and AA sepa- 
rated by four critical points —{2N — 1)02^2/^, —CiPi — 
C2/32, ci/3i-f C2/32 and (27V-l)c2;32/(A^ + l) corresponding 
to C12/3 -4, -3, 3, and 4 in Fig. [Tb. 
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In the region C12P G {-{2N -I)c2l32/N, -ci^i -C2/32), 
the MM_ phase takes the same form as eq. (|15p . with 
extreme states represented by 



\F,,F2,F, F) ^ 0^%%"^;/^ \N, N) \F2,F2) 



(18) 



or Zi/2(a|)^(&|)^^(Bt)(^-^^=^)/2 |0). In MM_ phase the 
atoms in species 1 are totally polarized in one direction 
and form a "steady magnetic field" (black solid line), 
and those in species 2 are partially polarized in the same 
direction. Increasing the strength of coupling interaction 
(|ci2/?|) breaks singlet pairs in species 2 one by one, and 
results in the increase of the total spin. 

In the region C12/3 G (— ci/3i— C2/32, ci/3i+C2/32), the two 
species are essentially independent for a weak interspecies 
spin-exchange interaction. The PP phase is a total spin 
singlet described by the direct product of the well known 
polar ground state Z"^ (it)A^/2(5t)A'/2 |o) giving 

rise to (Pfj = 0, ^F^^ = 0, and ^Fi • F2) = 0. 

In the region of ci2/3 G (ci^i +C2/32, (2iV- l)c2/32/(A^+ 
1)), however, the MM+ phase favors that the atoms in 
species 2 are polarizing to the opposite direction of those 
in species 1 and the total spin gradually decreases. This 
situation is much more complicated because all states 
that satisfy the condition mi +m2 = N — F2 are involved 
and the summation index runs over all possible Clebsch- 
Gordon coefficients, giving rise to the following extreme 
state 
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C. The case ci/3i < 0, C2/?2 > 



that in FF phase all atoms are polarized in the same di- 
rection and the total spin reaches its maximum value. 
The ground state is highly degenerate with degeneracy 
2F + 1. Take the state ^2, -F, = 0) as an example, 
the atomic populations are 
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which reduces to (iV/4, iV/2, iV/4) for each species for 
large N . The number fluctuations 
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are in order of vTV for large iV. On the other hand, 
in A A phase, atoms in species 1 and 2 are fully polar- 
ized, however, in opposite directions and the total spin 
vanishes [271. The singlet ground state is a fragmented 
[27I and the single particle density matrix 



condensate 



is diagonal with atoms equally populated (jIq^^j^ = N/S 
or (A^/3, N/3, N/3) for each species. The fluctuations 
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are anomalously large (in order of N) thus the frag- 
mented condensate is fragile. 



In the case of a mixture of a ferromagnetic and a po- 
lar condensate, four possible phases FF, MM_, MM+, 
and AA are separated by three critical points —{2N — 
I)c2l32/N,0 and (2iV - l)c2/32/(A^ + 1) corresponding to 
C12/? = -4, 0, 4 in Fig. 1 (c). 

It has been shown that the ground state of polar atoms 
(c2/32 > 0) only is a fragmented condensate with anoma- 
lously large number fluctuations and fragile stability 8] , 
which can be described as (Bt)^/2 |o). The B'< and B 
are invariant under any rotation of the system, and com- 
mute with F2 and i^2z- However, for a ferromagnetic 
condensate (ci/3i < 0) the ground state favors that all 
atoms are aligned in the same direction (i.e., {a\)^ |0)) 
and much more stable. So when these two kinds of atoms 
are mixed together, the polar atoms are more easier to 
be influenced, but their back action on to the stable fer- 
romagnetic atoms is negligible. This can be seen from 
the constant black solid line in FiglT](c). 

An interesting observation is that for large and neg- 
ative (positive) value of parameter C12P the system en- 
ters the same phase FF (AA), no matter how the atoms 
interact inside each species. Detailed calculation show 



IV. THE SITUATION OF 7 / 

The last term C127 in Hamiltonian ([5]) is related to the 
inter-species singlet paring. Although [F2^0|2Oi2] = 0, 
we notice that [F?, 812612] 0, and [Fi,e|26i2] 7^ 0. 
Thus in general they do not belong a set of commuta- 
tive operators and the system is not solvable. In the 
special case of ci/3i=C2/32=^Ci2/3, we found that the spin- 
dependent Hamiltonian reduces to a sum of commutative 
operators. The eigenstates can be constructed by sev- 
eral building blocks [13, [3 [IE Ull which can be found 
via generating function method. To see more clearly the 
role played by the 7 term, we focus on a special case in 
which C127 is much larger than the other parts of the 
Hamiltonian and ci/3i < 0, C2/32 > 0. We know that the 
7 term encourages pairing two different types of atoms 
into singlets when 7 < 0. In this case the interspecies 
singlet-pairing interaction dominates the system and the 

total spin vanishes, namely (^^^^ — ^(Fi +1^2)2^ = 0. 
There are two typical cases obeying the above condi- 
tion. In one case we have (f2\ = 0, (F2) = 0, and 
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FIG. 2: (Color online) The dependence of ground-state order 
parameters on ci2/3 at fixed values of ci/?i = —1, c^fi^ = 2, 
and C127 = —20 (in the unit of |ci/3i|). Black solid lines, 
red dashed lines, blue dot-dashed lines and green dotted 

lines denote respectively the order parameters (Fj 



Fi -Fa), and l?>. 



Fi • = 0, which means the singlet formation occm^s 

inside each species, and the atoms in the same species 
are all paired with no net spin left. In the other case we 

have -2 ^Fi • t-^ = {^Ff ^ + (tf^ ^ 0, which means 

some intra-species pairs are broken, meanwhile singlet 
formation occurs between different species. Under the 
condition Ni = N2 = 100, numerical results show that 
the above two cases indeed exist and well determined by 
the parameter ci2/3. 

We consider the direct product of the Fock states of 



t-u + • (1) (1) (1) 

the two species n\ , rig , n_{ 

may be equivalently defined as 
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Here mi_2 are the corresponding magnetization specified 

as TOi^2 = n^i''^^ — n^li"^ and m = mi + m2 is the total 
magnetization. For simplification, we restrict ourselves 
into the subspace that the total magnetization is con- 
served m = 0, in which case all states are non-degenerate. 
The Hamiltonian (|8]) is then represented in a sparse ma- 
trix and with the approach of exact diagonalization we 
numerically get the ground state of the system, on which 
the order parameters are calculated. Fig. 2 shows the 
dependence of the four order parameters on ci2/3 at fixed 
values of ci^i = —1, 02^2 = 2, and C127 = —20. We 
find that in the region of ci2/3 < 0, (F^V (Fj), and 



Fi -F; 



are all approximative ly equal to zero, corre- 
sponding to the first case of vanishing total spin. How- 
ever, for positive ci2/3, atoms in each species begin to 




FIG. 3: (Color online) The dependence of the four order pa- 
rameters on the parameter C12/? for different inter-species par- 
ing interaction parameter C127 = 0, —5, —10, —15, —20. 



polarize and are fully polarized for ci2/3 > 2, with the 
negative value of ■ ¥2'^ showing that they are polar- 
ized to the opposite directions. During the whole process 
the three order parameters always obey the condition 

-2(Fi- F2) = ( F?\ + ( fi). The distinction in the 



envelope of (fij and {flj ™ the region < ci2/3 < 2 
lies in the choice of ci/3i = —1, C2/?2 = 2. 

The fourth order parameter ^B|2012^ /3 refiects the 
feature of the total spin F of the system. It takes a 
constant value, which is actually the maximum value, 
in entire region of C12P in Figure 2 (see the green line), 
while the total spin reaches its minimum value 0. In 
order to see more clearly the competition between the 
interspecies coupling and singlet pairing interactions, in 
Figure 3 we illustrate four order parameters for differ- 
ent values of C127. An obvious variation can be easily 
seen from, for example. Figure 3d, where the black solid 
line (7 = 0) gradually changes to the purple dashed line 
(7 — —20). We find that when 7 =0 the order parameters 

1^X2^12^ /3 becomes zero in the FF phase implying that 
the total spin F gets to its maximum (Fig. 3a), while in 
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the AA phase it equals to a constant value and the total 
spin F vanishes. In between the two limits "0" and "con- 
stant", the intermediate value of the order parameters 
^€)|20i2^ /3 indicates that the system is a mixture with 
both singlet pairs and nonzero net magnetization. This 
feature can be well understood for the special case when 
7 = and /3 = 0, marked as a small dip in the black solid 

line in Fig. 3a. We find that (el20i2)/3 = 10000/9, 

which can be obtained analytically. In fact, the state 
at this point can be simply expressed as a direct prod- 
uct of two well-known states (ferromagnetic and polar), 
(a|)^ |0) (g) (i3t)JV/2 |Q^^ gj^j ^i^j^^ ^Yns state 



©12012) = 



- (a|ai6'_i&-i 



3 ' 



which agrees with the numerical result. We notice that 
the system is a mixture with N /2 polar pairs and net 
magnetization F = Fi = N. 

In our system there exist three states, on which the 
average of total spin amounts to zero. Although we 
have F2(^t)A^/2(5t)JV/2 |o) = o, F2 \n,N, 0,0) = 0, and 

^^(©12)^ |0) = 0, the difference between these states can 



be easily seen from an example Ni — N2 = N 
find that 



2. We 
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+C^2°:i°2,-i(A-)^-^(F2-)^+^(aI)^(6l 

+c'4,,,^,{F^^rHF2-rHa\rib\ 
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From the relation between these three states (|24l) , we see 
that the AA phase of our system includes at lest two 
pairing mechanism, i.e. and A^B'' . As total spin F 
vanishes, the number distributions of these three states 



,0) 



I — yt_i^ — N/3, but the number 

fluctuation on these states are quite different. For the 
state (it)JV/2(st)A'/2 |o) [i], the results are 
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On the state jA'', iV, 0,0) [23], we have obtained in eq. 
([22l) . in the case of large N, satisfying 



An 



2 ( An, 



An'/i 
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while for the state (QI^)'^ |0) , we find that the number 
fluctuations are equally distributed [131, i.e. 



An^^^) = (An^^'^) = (AnL^]) 
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(27) 




10000 



E 
to 
i_ 

to 

Q. 
i_ 
0) 
T3 



— <f;> 

- -<F> 

<F^.F2> 
<e-e>/3 



(b) 





C12P 



FIG. 4: (Color online) (a) The dependence of the fourth order 
parameter on the parameter C12P for different inter-species 
paring interaction parameter C127 — 0,1,5,10,20. (b) The 
dependence of ground-state order parameters on ci2/3 at fixed 
values of ci/3i = —1, C2P2 = 2, and C127 = 20. 

Finally we discuss briefly the situation when C127 > 0, 
in which case the ground state favors minimizing the 7 
term. We illustrate the competition between the inter- 
species coupling and singlet pairing interactions numer- 
ically in Fig. 4a. If C127 is far more larger than any 

other parameters, the order parameter ^0|20i2^ /3 will 

vanish (see the dashed line in Figure 4a). All order pa- 
rameters take the value or ~ N'^ when ci27 ^ 1, with 
the boundary determined by the corresponding ampli- 
tudes of spin-coupling strengths as shown in Fig. 4b for 
C127 = 20. 
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V. CONCLUSION 

In summary we studied the interspecies singlet pairing 
in the ground state of a binary mixture of spin-1 conden- 
sates in the absence of a magnetic field. In the case of 
C127 = 0, the exact quantum states can be constructed 
from angular momentum theory for the mixture of two 
ferromagnetic, two polar, and ferromagnetic-polar con- 
densates. The ground state is classified into five types 
according to the inter-species coupling parameter ci2/3. 
By means of the full quantum approach of exact diago- 
nalization, more general case of 7 ^ is considered. We 



illustrate the competition between the two interspecies 
interaction Ci2/3 and C127, and find that if C127 <t; — 1 
the ground state is a singlet of the total spin. There, 
however, exist different types of singlet formations deter- 
mined by C12/3. 

This work is supported by the NSF of China under 
Grant Nos. 10774095 and 11074153, the National Basic 
Research Program of China (973 Program) under Grant 
Nos. 2010CB923103, 2011CB921601, the NSF of Shanxi 
Province, Shanxi Scholarship Council of China, and the 
Program for New Century Excellent Talents in University 
(NCET). 



[1] D. M. Stamper-Kurn, M. R. Andrews, A. P. Chikkatur, [21] 
S. Inouye, H.-J. Miesner, J. Stenger, and W. Ketterle, 
Phys. Rev. Lett. 80, 2027 (1998). [22] 
M. Ueda and Y. Kawaguchi, arXiv: 1001.2072 
H.-J. Miesner, D. M. Stamper-Kurn, J. Stenger, S. In- 
ouye, A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. [23] 
82, 2228 (1999). 

T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 [24] 
(1998). 

C. K. Law, H. Pu, and N. P. Bigelow, Phys. Rev. Lett. 
81, 5257 (1998). 

H. Pu, C. K. Law, S. Raghavan, J. H. Eberly and N. P. 
Bigelow, Phys. Rev. A 60, 1463 (1999). 
W. Zhang, D. L. Zhou, M-S. Chang, M. S. Chapman, 
and L. You, Phys. Rev. A 72, 013602 (2005) 
T.-L. Ho and S.-K. Yip, Phys. Rev. Lett. 84, 4031 (2000). [25 
T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

E. J. Mueller, T.-L. Ho, M. Ueda, and G. Baym, Phys. [26 
Rev. A 74, 033612 (2006). 

J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Mies- [27 
ner, A. P. Chikkatur, W. Ketterle, Nature (London) 396, 
345 (1998). [28 
M.-S. Chang, C. D. Hamley, M. D. Barrett, J. A. Sauer, 
K.M. Fortier, W. Zhang, L. You, and M. S. Chapman, [29 
Phys. Rev. Lett. 92, 140403 (2004). 

M. Koashi and M. Ueda, Phys. Rev. Lett. 84, 1066 [30 
(2000). 

M. Ueda and M. Koashi, Phys. Rev. A 65, 063602 (2002). [31 
C. V. Ciobanu, S.-K. Yip, and T.-L. Ho, Phys. Rev. A 
61, 033607 (2000). [32 
T.-L. Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 3276 
(1996). [33 
H. Pu and N. P. Bigelow, Phys. Rev. Lett. 80, 1130 
(1998); ibid, 80, 1134 (1998) 

B. D. Esry, C.H. Greene, J. P. Burke, Jr., and J. L. Bohn, [34 
Phys. Rev. Lett. 78, 3594 (1997). 

E. Timmermans, Phys. Rev. Lett. 81, 5718 (1998). [35 

C. J. Myatt, E. A. Burt, R. W. Christ, E. A. Cornell, 
and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997). 



G. Modugno, M. Modugno, F. Riboli, G. Roati, and M. 
Inguscio, Phys. Rev. Lett. 89, 190404 (2002). 
G. Thalhammer, G. Barontini, L. De Sarlo, J. Catani, F. 
Minardi, and M. Inguscio, Phys. Rev. Lett. 100, 210402 

(2008) . 

S. B. Papp, J. M. Pino, and C. E. Wieman, Phys. Rev. 
Lett. 101, 040402 (2008). 

P. O. Fedichev, Yu. Kagan, G. V. Shlyapnikov, and J. 
T. M. Wahaven, Phys. Rev. Lett. 77, 2913 (1996); F. 
K. Fatemi, K. M. Jones, and P. D. Lett, ibid. 85, 4462 
(2000); M. Theis, G. Thalhammer, K. Winkler, M. Hell- 
wig, G. Ruff, R. Grimm, and J. H. Denschlag, ibid. 93, 
123001 (2004); C. D. Hamley E. M. Bookjans, G. Behin- 
Aein, P. Ahmadi, and M. S. Chapman, Phys. Rev. A 79, 
023401 (2009). 

Z. F. Xu, Y. Zhang, and L.You, Phys. Rev. A 79, 023613 

(2009) . 

Z. F. Xu, J. Zhang, Y. Zhang, and L. You, Phys. Rev. A 

81, 033603 (2010). 

J. Zhang, Z. F. Xu, L. You and Y. Zhang, Phys. Rev. A 

82, 013625 (2010). 

Z. F. Xu, J. W. Mei, R. Lii, and L.You, Phys. Rev. A 
82, 053626 (2010) . 

M. Luo, Z. Li, and C. Bao, Phys. Rev. A 75, 043609 
(2007). 

W. Zhang, B. Sun, M. S. Chapman, and L. You, Phys. 
Rev. A 81, 033602 (2010). 

A. C. Maan, E. Tiesinga, H. T. C. Stoof, and B. J. Ver- 
haar, Physica. B. 165, 17 (1990). 

S. B. Weiss, M. Bhattacharya, and N. P. Bigelow, Phys. 
Rev. A 68, 042708 (2003); ibid. 69, 049903(E) (2004). 
A. Pashov, O. Docenko, M. Tamanis, R. Ferber, H. 
Knoeckel, and E. Tiemann, Phys. Rev. A 72, 062505 
(2005). 

S. Yi, 6. E. Mustecaplioglu, C. P. Sun, and L. You, Phys. 
Rev. A 66, 011601(R) (2002). 

T.-L. Ho and L. Yin, Phys. Rev. Lett. 84, 2302 (2000). 



